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£N| , We argue that the < A OPE > condensate found in the Landau gauge on lattices, when an Operator 

Product Expansion of Green functions is performed, might be explained by instantons. We use 
cooling to estimate the instanton contribution and extrapolate back the result to the thermalised 
configuration. The resulting < Af nst > is similar to < Aq PE >. 
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■ I- INTRODUCTION 

Lattice calculations of the gluon propagator and three-point Green function in the Landau gauge indicate that the 
^v^j [ expected perturbative behaviour at large momentum p has to be corrected by a 0{l/p 2 ) contribution sizeable up to 
10 GeV [0-|5|. An understanding of this contribution as the effect of an A 2 = A®A£ condensate (in the Landau gauge 
A 2 is the only mass dimension-two operator liable to have a v.e.v.) has been gained by verifying that two independent 
Green functions could be described by the perturbative contribution corrected by the effect of one common value of 
< Aq PE >, as expected from OPE. The physical origin of this condensate is an important question, possibly involving 
the non-trivial topology of the QCD vacuum. In particular, instantons provide an interesting insight into a wide range 
of low energy QCD properties ( M and refs therein). They have been put into evidence on the lattice using different 
cooling procedures. In this letter we claim that instantons provide for < A 2 > a value close to what is needed for the 
OPE fit to Green functions. 

We propose a method to identify instantons from the cooled gauge configuration, count them and measure their radii; 
we also check that these results are compatible with the instanton number deduced from the two-point correlation 
function of an instanton. We then estimate < ^4 2 nst >, the contribution of the instantons to < A 2 > in cooled 
configurations, extrapolate back to the thermalised configurations (zero cooling sweeps) and finally, compare the 
outcome with the OPE estimate. 



II. COOLING AND INSTANTON COUNTING BY SHAPE RECOGNITION 



A. cooling 

In order to study the influence of the underlying classical properties of a given lattice configuration, the first step 
will be to isolate these structures from UV modes. The method we use is due to Teper Hj; it consists in replacing 
each link by a unitary matrix proportional to the staple. A cooling sweep is performed after replacing all the links of 
the lattice. This procedure introduces largely discussed biases, such as UV instanton disappearance and instanton- 
anti-instanton pair annihilations, that increase with the number of cooling sweeps; alternative cooling methods have 
been proposed (see for example || and refs. therein) to cure these diseases. We will try to reduce them by identifying 
the instantons after a few cooling sweeps and extrapolating back to the thermalised gauge configuration. 
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B. Instantons 



Instantons (anti-instantons) are classical solutions |j of the equations of motion. We work in the Landau gauge 
which is defined on the lattice by minimizing A^(x) 2 . For an instanton solution this prescription leads to the 
singular Landau gauge, where the gauge field is 



A (I)a = 27 ? ° y 3:„p 2 

" x 2 (x 2 +p 2 ) ' 



p being the instanton radius. The instanton has been chosen centered at the origin with a conventional color orien- 
tation. These solutions have Q = ±1 topological charge, 

Q 55 32^ / ** F ir%» = / d A xF^F; v , (2) 

where F^ v = \tp,vpaF pa . From eq. (|l|) the topological charge density is: 



2 x 4 



On the lattice, the topological charge density will be computed as 

±4 



0) = ^^^2^^ P aTr[U^(x)U pa (x)} , (4) 



Hatt 

±1 

with IL^^a;) = L^(s)^(x + /ia)L^(a; + /ia + ^a)C/t(x + ^a), and ? Al „ pcr the antisymmetric tensor, with an extra minus 
sign for each negative index. 

C. identification of instantons 

A common belief is that an instanton liquid gives a fair description of important features of the QCD vacuum. 
Along this line we will try a desciption of our cooled gauge configuration as an ensemble of non-interacting instantons 
with random positions and color orientations. We hence also neglect the interaction-induced instanton deformations 
and correlations. Although the instanton ensemble for the QCD vacuum cannot be considered as a dilute gas — p~2]| , 
this crude assumption allows a qualitatively reasonable picture, especially near the instanton center. 

Many enlightening works have studied the instanton properties from lattice gauge configurations^, among which 
[fl3| |l6| . As for us, we start by searching regions where the topological charge density looks like that of Eq.(||). 
Starting from each local maximum or minimum of Qiatt(^) we integrate over all neighbouring points with |Qi a tt(a:)| > 
<x\Qia,tt(%max)\) for different values of a ranging from 0.8 to 0.4. A local extremum is accepted to be an (anti-)instanton 
if the ratio e between the lattice integral and its theoretical counterpart, Q p (x), 

1 - 3a 1 / 2 + 2a 3 / 4 r 1 / d 4 x Q utt (x) (5) 

7 Jx/ |Q(x)|>q|Q(0)| 

shows a plateau when a is varied. Indeed for a theoretical (anti-)instanton e = 1 (e = — 1) for any a € [0,1]. 
As a cross-check of self-duality, this instanton shape recognition (ISR) procedure is applied on the lattice to both 
expressions for Q introduced in eq. (^|). 



X A detailed comparison of our results to their's will appear in a forthcoming extensive work. 
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D. Instanton numbers and radii 



The ISR method resolves the semiclassical structure on the lattice with only ~ 5 cooling sweeps. This early 
recognition reduces the possible cooling-induced bias. 

Whenever this method succeeds, we measure the radius of the accepted instanton in two separate ways: from the 
maximum value of Q, written QJ^j*, at the center of the instanton, and from the number N a of lattice points in the 
volume f , ,„,„,, d A x. We find: 

Jx/ |Q(a;)|>a|Q(0)| 



'latt 

These two measures of the radius agree within 10%. 



Va-V4 _ i 



P/ a = \ onmax = / i/„ : V — > ( 6 ) 



III. GLUON PROPAGATOR 
A. The instanton gauge-field correlation function 

The classical gauge-field two-point correlation function verifies, for any position and color orientation, 

lUKl = ^£(4' )a (^ )a (-*)) -G (2) ( fc 2 ) (V-^f) , (7) 

a a ^ ^ 

where V is the volume in the cuclidean four-dimensional space and A^ )a (k) is the Fourier Transform of Eq. ([[]). The 
resulting scalar form factor is: 



,( 2)/ ,,. _ 327T 4 ^ (hp) 2 



K% being a Bessel function JlTj . Eq. (||) equally applies to instantons and anti-instantonsP]. 

In a perfect gas approximation for an ensemble of ni (jia) (anti-)instantons of radius p, the classical gauge-field 
correlation function is simply given by Eq.(|^) times the number of instantons and anti- instantons, nj + tia- This 
correlation function is the contribution of the background field to the gluon propagator. We expect this formula to 
describe the behaviour of the lattice gluon propagator once the effect of quantum UV fluctuations is removed by the 
cooling procedure Q. The effect of instanton interactions is known ||l0| , |lT|| to modify the instanton shape far from 
its center, in the IR region. But the large k 2 behaviour should be appropriately given by Eq.(||) i.e. oc 1/fc 6 . This is 
shown in Fig.|l|.a for one generic lattice gauge field configuration. 

The theoretical lines in that plot ^] are generated by Eq.(||) using the average radius p and ni + ua computed 
from the ISR method. Note that the matching improves with the number of cooling sweeps. This agrees with the 
expectation that decreasing the instanton density reduces the instanton deformation and that quantum fluctuations 
are damped by cooling. Reversely, if we know the average radius from the ISR method, we can compute ni + ha from 
the fit to the measured propagator. 

B. The hard gluon propagator 

Let us now consider a hard gluon of momentum propagating in an instanton gas background. The gluon 
interacts with the instanton gauge field. This can be computed with Feynman graphs and it is easy to see that when 
the instanton modes verifies fc M ~ 1/ p <§C p^, the dominant contribution is an 0(l/p 2 ) correction to the perturbative 

gluon. This correction is equal to the standard OPE Wilson coefficient [|[D times < Af nst (0) >= 1/V J d 4 x (A^ 1 ^) 2 (x). 
We will now proceed to estimate this instanton-induced condensate. 



2 A similar analysis is being parallelly performed, in other context, by Broniowski and Dorokhov |T^] 
3 We multiply by k 2 to compute a dimensionless object and perform the matching. 
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IV. < A 2 > CONDENSATE 



A. < A > in instantons 

From Eq.(|l|) we get 

< AIM >- / ^EW" - 12^^^ . (9) 

where p is the average instanton radius in the considered cooled configuration and n c is the number of cooling sweeps. 

We use an ensemble of 10 independent gauge configurations^ at (3 = 6.0 on a 24 4 lattice. Each configuration has been 
cooled and after 5,7,10,15,30 and 100 cooling sweeps transformed into the Landau gauge. Using the ISR method on 
each gauge configuration, we obtain the results of Tab.|. In this table we also present the number of (anti-)instantons 
and the corresponding value for < A 2 (n c ) > computed by a Correlation Function Fit (CFF) i.e. a fit of the lattice 
propagators to the instanton correlation function, Eq.(||). The CFF method is expected to be affected differently 
from the ISR method by systematic uncertainties: instanton interactions, deformations and quantum fluctuations (as 
we can see in Fig.[l].a, at low momentum), and we therefore consider as quite encouraging the qualitative agreement 
- becoming quantitative at large n c - between ISR and CFF results. We have, for simplicity, translated our lattice 
results into physical units using, for all values of n c , the n c = inverse lattice spacing, a~ 1 (n c = 0) = 1.996 GeV (at 
P = 6.0). This simple recipe overlooks the effect of cooling on the lattice spacing (see ref. |l5| and refs. therein) but 
this simplification becomes harmless after extrapolating back our results to n c — 0. 



B. < A 2 nBt (n c ) > at zero cooling. 



The instanton number depends on the number of cooling sweeps. This result may imply that the cooling procedure 
destroys not only quantum UV fluctuations but something else from the semiclassical background of gauge fields. 
To lessen this problem we take advantage of the early recognition of the instanton content in a gauge configuration 
ensured by the ISR method and perform an extrapolation JL9| to n c = of the ISR results for < Af nBt (n c ) > in the 
table. We then obtain (see fig.l.b): 

< Al lst (n c = 0) >= 1.76(23) GeV 2 . (10) 

We have used a form a/(b + n c ) to fit and extrapolate. We have also varied a little this functional form to check 
the stability of the extrapolation. We take this result as indicative of the non-perturbative instanton contribution to 
the < A 2 > condensate. If we applied other lattice estimates of instanton gas parameters taken from the available 
literature to Eq. (|9]), the value of < A^ > would range [] from 1 to 2 GeV 2 . On the other hand, parameters from 
instanton liquid based phenomenology Q yield estimates of the order of 0.5 GeV 2 . As the quoted error in Eq.fjll]) is 
only statistical, this last range somehow estimates a certain systematic uncertainty. 





P (fm) 


ni 


+ riA 


< A 2 > 


(GeV 2 ) 


ISR 


ISR 


CFF 


ISR 


CFF 


5 


0.329(2) 


87(2) 


93(10) 


1.38(3) 


0.9(8) 


7 


0.361(2) 


74(2) 


59(1) 


1.42(5) 


1.12(2) 


10 


0.394(4) 


60(1) 


38(1) 


1.36(4) 


0.86(2) 


15 


0.417(5) 


43(1) 


28(1) 


1.11(4) 


0.72(2) 


30 


0.452(9) 


26(1) 


19(1) 


0.80(6) 


0.57(3) 


100 


0.53(1) 


10(1) 


9(1) 


0.43(3) 


0.37(3) 



TABLE I. Estimates of < A 2 nBt (n c ) >. 



4 Considering the present size of our systematic uncertainties we did not consider it worthwhile to increase further the statistics. 

5 We use the parameters obtained in for simulations on different lattices and f3's with a cooling improved to let scale 
invariant instantons solutions exist for large enough instanton sizes. We only quote anyway the results where the packing rate, 
n 2 /2(p/L) 4 (nj + ua), is as much as 1, since our method to estimate < A? nst (n c ) > assumes limited overlap between instantons. 
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FIG. 1. In fig. (a) (left) we present the lattice gluon propagators after several number of cooling sweeps (points) and the 
corresponding theoretical instanton gauge- field correlation functions (lines) in the perfect instanton gas approximation, Eq.(|^), 
plotted as a function of n 2 (n M = L/(2ir) L being the lattice length). In fig. (b) (right), we show the extrapolation at zero 
cooling of < A 2 (n c ) > in physical units (for three different trial functions). 



C. Comparison with < A > from OPE. 



Our instanton estimate of (A 2 ) is a semiclassical one, deprived of the necessary UV fluctuations, and therefore not 
directly comparable with ||] (AQ PE )(10GeV) = 2.4(5) GeV 2 . There is of course no exact recipe to compare both 
estimates, since the separation between the semiclassical non pcrturbative domain and the perturbative one cannot 
be exact. We may appeal to the fact that at the renormalisation point /i, the radiative corrections are minimised; 
therefore a semiclassical estimate must best correspond to (Aq PF ) at some reasonable /i. In the example of <j> 2 vacuum 
expectation in the spontaneously broken r/> 4 model given in pi[ , one finds indeed that it equals the classical estimate 
for fi around the spontaneously generated mass. In our problem, one could guess that the corresponding scale should 
typically be around l/p — 0.7 GeV, or some gluonic mass, a very low scale anyway. We cannot run (AQ PE )(10GeV) 
down to such a low scale 101, 



A 2 - A 2 



1+35 fcil/k /L_ 
44 no Amom 



fill 



we therefore stop arbitrarily /i ~ 2.6 GeV, where the OPE corrected perturbative running of the Green functions fails 
to correctly describe their behaviour. This scale of 2.6 GeV turns out to be of the same order [|| as the critical mass 
p2| for the gluon propagator. At this scale, we obtain: 



< ^qpe(2-6 GeV) >= 1.4(3)(3) GeV 2 , 



(12) 



where the first quoted error just propagates the uncertainty from the OPE determination of < A 2 > and the second 
one takes into account, in the way proposed in ref. Q, higher orders^ in a s for running. 

V. DISCUSSION AND CONCLUSIONS 

We are aware that our method of comparison of < ^4 2 nst (n c = 0) > and < Aq PE > suffers from a lot of arbitrari- 
ness and approximations (such as the perfect gas approximation, possible errors in the instanton identification, the 
uncertainty in the extrapolation to zero cooling sweeps, etc.). We have taken care to crosscheck our estimates by 
comparing different methods at each step of the computation, in particular the ISR and CFF (see Fig.|[a and Tab. I). 



3 We re- write Eq.Q in terms of the coupling constant renormalised in other schemes, such as MS. 
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A comparison with direct "measurements" of the (A 2 ) condensate from cooled lattice configurations could be thought 
as an additional crosscheck. Qualitative agreement is found for a large enough number of cooling sweeps, but this 
agreement is manifestly destroyed by UV fluctuations already for n c ~ 30. Of course, by using ISR and instanton gas 
approximation we sharply separate UV fluctuations from the semiclassical background. All these imprecisions seem 
anyway inherent to the subject. 

With this in mind, we nevertheless take the fair agreement between Eqs. ( |l0| ) and ([l2]) as a convincing indication 
that the A 2 condensate receives a significant instantonic contribution. In other words, the instanton liquid picture 
might yield the explanation for the 1/p 2 corrections to the perturbative behaviour of Green functions computed with 
thermalised configurations on the lattice. 
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